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We study a nonsingular bounce inflation model, which can drive the early universe from a contract¬ 
ing phase, bounce into an ordinary inflationary phase, followed by the reheating process. Besides the 
bounce that avoided the Big-Bang singularity which appears in the standard cosmological scenario, 
we make use of the Horndesky theory and design the kinetic and potential forms of the lagrangian, 
so that neither of the two big problems in bouncing cosmology, namely the ghost and the anisotropy 
problems, will appear. The cosmological perturbations can be generated either in the contracting 
phase or in the inflationary phase, where in the latter the power spectrum will be scale-invariant 
and fit the observational data, while in the former the perturbations will have nontrivial features 
that will be tested by the large scale structure experiments. We also fit our model to the CMB TT 
power spectrum. 


I. INTRODUCTION 

The standard model of cosmology regards inflation [H-Q as an important period in the early universe. As a super¬ 
fast expansion after the Big Bang, the inflation can solve a series of cosmological problems such as horizon, flatness, 
monopole and so on, as well as give rise to scale-invariant scalar perturbations that can fit with the data. However, 
traditional inflation scenario still needs to be improved, because of the so-called “singularity problem” which was 
proposed by Hawking et al. in their early work of singularity theorem [3, § . According to their proof to this theorem, 
if we track back inflation to its very beginning, we will generally meet the singularity of the early universe (the Big- 
Bang singularity). At the singularity, everything blows up and one can not get control of the universe under classical 
description. Since the singularity occurs before the onset of inflation, it can hardly be solved within inflation scenario 
itself. This motivates us to find alternative theories in pre-inflation era. 

Phenomenologically, there might be quite a few evolutions that can be set in front of inflation in order to get 
rid of this problem. As an example, the universe may undergo a contracting phase where the scale factor a(t) 
shrinks initially and then, by some mechanism, “bounces” into an expanding one (9| . The whole process can be done 
non-singularly if at the bouncing point aB{t) yf 0 M- The bouncing scenario has many interesting properties, for 
example, the Big-Bang puzzles such as horizon problem and flatness problem can be solved even in contracting phase, 
and scale-invariant primordial perturbations can be generated, etc [ll|. Moreover, such non-singular scenario can also 
be non-trivially extended to the cyclic universe [IJ, [3 • 

In bouncing inflation scenario however, there will be two more latent problems. One of them is the so-called 
“ghost instability”. This problem claims that field theory models that violates the Null Energy Condition (NEC) will 
generally give rise to “ghost” degrees of freedom, which causes the instability. The problems is originated from the 
“Ostradski problem” in field theory, and the first cosmological application of this problem was on dark energy 
models (see also |l6j| . which discussed from quantum aspect). Since a nonsingular bounce violates NEC, normally 
the ghost will appear. However, recently an interesting “Galileon” theories [13) [3 been shown to be able to avoid 
such ghost degree of freedom while violating NEC (see pioneering work on such kind of models in 0)- The reason 
is that although Galileons have higher derivative operators that can violate NEG, due to the delicate design of the 
Lagrangian, there will only be one dynamical degree of freedom, which can be made normal. The other degrees of 
freedom is non-dynamical, thus will not lead to any instabilities. Therefore, in this paper we will make use of Galileon 
theory to build our model. Note that pure Galileon bounce models has also been proposed in the literatures 0 , 111 , 
which can avoid ghost instability on NEC violation. 

The other problem is known as “anisotropy problem” (H , which will impose additional constraint on the evolution 
of the universe before the bounce. An exact isotropic universe at the initial time needs somehow fine-tuning, so in 
realistic models of the early universe, certain amount of anisotropies will exist. In pure inflation scenario, this will 
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not be a problem because the anisotropies will decay in expanding universe, and will eventually be diluted away by 
inflation to get an isotropic universe. In bounce cosmology however, the anisotropies will grow in contracting phase, 
and if the growth of the background energy density does not exceed that of anisotropies, the latter will dominate the 
universe to make it a totally anisotropic one, and probably no bounce will occur and the universe will never enter into 
an expanding phase, which shows contradiction to our observations. By calculation, we can know that the growth of 
anisotropies scales as a“®, and to have growth of background faster than this speed one must have the equation of 
state (EoS) larger than unity, which gives constraint to the model building. In fact, a large EoS can be obtained by 
requiring a negative potential. Though currently such a potential is only for phenomenogically use, it can actually be 
reduced from more fundamental theories such as Ekpyrotic theories. 

The phenomenology of bounce inflation scenario was first studied in [2^ , where it was shown that such a scenario can 
not only obtain scale-invariant scalar perturbations that can fit with the data, the perturbations generated before the 
bounce can also give rise to tilted spectrum, which can explain the suppression of CMB TT spectrum at large scales. 
In the second paper of Ref. 0 , a bounce inflation model was realized, but the ghost and anisotropy problems has 
not been addressed. In [2^, a sketch plot of healthy bounce inflation model has been drawn with the ways of avoiding 
those problems mentioned above, and in [2^, the Ist-year PLANCK data has been used to constrain this scenario. In 
(also see [ 13 ), it was found that the bounce inflation may generate a large circularly polarized gravitational wave 
on large scale, which will lead to the enhancement of TB and EB-mode correlations on corresponding scale, see [l^l 
for recent study on large scale power deficits on the CMB fluctuations along with various pre-inflationary models. 

In this paper, we aimed at building a realistic bounce inflation model, and fit with the newest PLANCK data, 
hollowing 2^ and also the above arguments, we make use of the Galileon theory the avoid ghosts, and introduce 
negative potential to keep anisotropies away during the evolution. Moreover, another technical point is how to unify 
a negative potential with a positive potential which is used to drive a period of inflation, and in the following we 
will show that, with the help of some “shape functions”, one can naturally match the evolutions of contracting and 
expanding phases together. In the case of model building the shape functions are shown to be very useful, although 
its fundamental origination needs to be more explored. 

The rest of the paper will be organized as follows: in Sec. II we set up our model, and show by numerics that it 
can give us non-singular bounce followed by a period of inflation, and the ghost and anisotropy problems will not 
appear. In Sec. Ill we calculate the perturbations generated in our model, and show that for perturbations of small 
wavelengths, nearly scale-invariance will be realized similar to inflation case, while for those of long wavelengths, there 
will be a blue-tilt which can explain the low-Z suppression in CMB TT spectrum, in Sec. IV we will fit our model 
with the newly released PLANCK data. Sec. V comes as a conclusion. 


II. THE GALILEON BOUNCING INFLATION MODEL 


As has been demonstrated in the introduction, in order to avoid ghost instability, it is useful to make use of Galileon 
theory to build up the bounce inflation model. The most generalized Lagrangian of Galileon is constructed as [l^l : 
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where X = —d^<j)d^(j)/2 is the canonical kinetic term of </>, □</> = R is the Ricci scalar and Gf^u = 

~ g^ivR/‘2‘ is the Einstein tensor. However, in order to make our model simpler, we will only use £2 and £3 
instead of all the £i’s. To be explicit, we take the Lagrangian to be: 

£ = k{(j))X + t{(j))X^ - V{y) - G(A, , (1) 


where the last term, which we called the “G-term”, contains the second derivative of the scalar field (j) which will be 
essential in violating NEG and trigger the bounce. Note that for the kinetic part k{(j))X + t{(j))X'^, when k{(j)) = 1, 
t{(p) = 0, it will reduce to trivial kinetic term X, and when k{(p) = —1, t{(p) = 1 it reduces to the form of “ghost 
condensate”. We will see in the following that kinetic part of this form is useful in matching bounce with inflation. 
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which was first proposed for ghost-free bounce by Cai et al [1^, and later connected with Supergravity in [30| . 
Following the Lagrangian o, one can get the equation of motion for (j) as: 

[k{4>) + Qt{(j>)X + &GxH^ + 6HGxx^ - 2(G^ + Gx^X)]^ + 37J[fc(<^) + 2t(<^)X - 2(G^ - Gx^X)]^ 

+ [2k4^) + 44(,^)X + 6Gx{H + 3H^) - 2G^^]X - k^^) - m(j))X^ + = 0 , (2) 

and the energy density and pressure are 

p = k{(j))X + 3t{cj))X^ + 3GxH<j)^ - 2G^X + V{<t)) , P = k{<j))X + t{cj))X^ - 2{G^ + Gx^)X - V{(j>) (3) 

respectively. 

We can furtherly determine the functions in Lagrangian o with the evolutions we want to get. First of all, an 
anisotropy-free contraction phase which requires the EoS larger than unity generally leads to a negative potential, the 
simplest Lagrangian of which is: 


^con ^ ^ 


( 4 ) 


where < 0. Meanwhile, the inflationary expansion phase after the bounce will have a positive flat potential: 

£“/=X-(5) 

where > 0, <C 1. Therefore, we will choose functions k{(j))^ G{X,(j)) and V{(j)) so that 

they can approach Lagrangians o and m in the limits of past and future, respectively. As a simple example, we 
choose forms of k{(j)), and G{X,(j)) as: 
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( 6 ) 


and the potential is: 

E(0) = [l-tanh(Ai^)]W°'‘(0) + [l + tanh(A2^)]E“^((/>) , = A4(1 - ^)2 , (7) 

Mp Mp 

with kg, t n, Ki. Ai, A 2 , Vqj c, A and v are all free parameters. Note that ((?!)) is inspired from the Ekpyrotic 
potential pll. [^ [^. and E*"'^((/>) is the well-known spontaneously symmetry breaking potential [s^ which fits well 
with the PLANCK data [s^. So the choice of this potential is not only simple, but may have a strong background of 
fundamental theories. 

From Eqs. ([6]) and ([2]) one can see that, in the {(j>/Mp) ^ 1 limit, we have 


fc(0) = l , t{cl>)=G{X,cj,) 
while in the {(j)/Mp) — 1 limit, we have: 

kicj,) = l, t(<^)=G(A,</.) 


0 , E((/>) = , 

(8) 

0 , Vi^) = , 

(9) 


which can give us (H]) and ([5]) respectively. Around (^ = 0 however, G{X, cj)) and the second term of fc(0) becomes 
large, so the linear term of X will flip its sign, while the A^-term and G-term will be dominate. Then for the solution 
where </> increases monotonically (which is a very natural solution), one will get the anisotropy-free contraction in the 
past and inflation in the future, between which the bounce is triggered around <j) = Q. 

Furthermore, we would like to emphasize that by having the shapes of k{(j))^ and G{X, <f), one can successfully 
control each term such that the higher derivative term only dominates for a while in order to trigger the bounce, 
while disappears before and after that, not disturbing the evolution before and after the bounce. In fact, if we do not 
control the higher derivative term, it will usually become more and more important after dominance, leading to a very 
fast evolution of the universe instead of a modest one, or even the “big-rip” singularity. Moreover, using tanh-like 
shape function, we can naturally connect and so that the field (j) can evolve naturally from a lower 

potential to a higher one, leading to inflation. Therefore we can see that, shape functions are actually very useful in 
constructing models. We also show the shape functions and potential explicitly by numerics in Fig. [1 and [H For 
pioneer works of using this shape functions to model building, see [ 2 ^ for ghost-free bounce and see 3^ for dark 
energy models. 
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FIG. 1: Plots of functions fc(<?f>), and X ^G{X,(j>) in Eq. ® w.r.t. the scalar, (j). Note that the last one is also function of 
(j> only. The parameters are chosen as ko — 0.6, ki = 15, to = 5, K 2 = 10, 7 = 1 x 10®. In such a choice, all the three functions 
have nontrivial value only around the bouncing point, which is useful to trigger the bounce. 


vi« 



FIG. 2: Plots of potential in Eq. ((2]) w.r.t. the scalar, 0. The parameters are chosen as Ai = A 2 = 10, Vo = 0.7Mp, 

c = Vw, A « 1.5 X 10~^Mp, V = lOMp. From the plot we can see that the shape of the potential shares the same shape of 
Ekpyrotic potential for negative (f) while that of symmetry breaking inflaton potential {(j)) for positive 4>. We will 

see below that 0 = 0 is almost the division of contracting and expanding phases. 


Now let’s go to a little bit more detail about the field evolution both before and after the bounce, which actually 
become very simple. In the contracting phase before the bounce when Eq. ([ 8 ]) applies, the equation of motion of (j) 
becomes: 


^ + = 0 . 


The energy density and pressure of ({> thus becomes 

pcon + l/“" = ~ ~ , 

while the EoS and the slow-roll parameter of the universe are 


pc 


Wc = 


~ 1 


(/-2 - 


3, , 6Vbe“'^/^r 


( 10 ) 


( 11 ) 


( 12 ) 


Since cj)^ > due to the requirement of positivity of the energy density, Wc is obviously larger than unity, 

and the growth of the energy density p oc will be larger than the anisotropy, avoiding the dominance of the 

latter. So the negative potential is rather useful in avoiding the cosmic anisotropy problem. 

Setting the initial values of (j) be negative, and once (j) approaches 0 with certain amount of velocity the higher 
derivative terms is turned on, which is responsible to trigger the bounce. After bounce however, as Eq. Q applies, 
these terms decay away again, and one will have: 

0 + 3iI0 + 2A44(^-l) = O , 

yZ yZ 


( 13 ) 
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which gives the solution of inflation. According to the potential, one can express the slow-roll parameter Ce as well as 
the e-folding number Nf. of the inflation as: 











(14) 


By requiring Ce(<()/) = 1 and Ne((j)) = 60, the field values at starting and ending points of inflation can be calculated 
as: 


Ipi Ri 1.9Mp , (j)f Ri 11.66Mp . 


(15) 


One can see again that the (j) has monotonically increasing evolution. The potential has valley with a minimum at 
V = 13Mp, and when (j) falls down the valley, it will become oscillating rapidly, which reheats our universe. The whole 
evolution of the (j), (j), H and w are sketched in Fig. [31 



FIG. 3: Top row: Plots of scalar (j) and its velocity cj) w.r.t. cosmic time t. During the evolution, (j) goes from negative to 
positive values, and as can be seen below, the switching point is nearly the bouncing point. The velocity peaks at the bouncing 
point, and decay towards both past and future directions. Bottom row: Plot of Hubble parameter H (left) and EoS w (middle 
and right). The bubble parameter transits from negative to positive values at the bouncing point, and will approach to certain 
constant positive value, indicating a de-Sitter (inflationary) expansion after the bounce. The middle plot is the EoS near the 
bounce region. From the plot we see that while before bounce it approach to value larger than unity, it goes down to minus 
infinity at the bounce point and approaches to —1 after the bounce. The right plot is the EoS during whole evolution, and since 
the inflationary period is very long, the features near the bounce is not viewable. Nevertheless, we can see that after inflation, 
EoS shows some oscillating behavior, referring to the reheating process. 


III. PERTURBATIONS 
A. scalar perturbations 

In this section, we analyze the perturbations that could be generated in this model. Following the general formulae 
for calculating perturbations of Galileon models in [H, [s^ , the 2-nd order perturbed lagrangian based on o turns 
out to be: 



drjdPxa 
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where 


Q 


2M^X 


[M^H - GxX(t>f 
{M^H - GxX^f 
2M^X 


[k{^) + 2t{4>)X + 2{Gx + GxxX)^ + AHGx^ - 

^^2 v2 

[fc{(^) + 6t{^)X + 6HiGx + GxxX)^ + ^ 


2G\X^ 

~Mf~ 

-^Q . 


(17) 

(18) 


Here C is the variable of adiabatic perturbation, while is the sound speed squared of representing how fast the 
perturbations can propagate. From the quantum field theories, the ghost-free condition of C is Q/c^ > 0. Moreover, 
the gradient stability requires 0 < < 1. From Eqs. di) and (fT8)l one can see that, for most of the evolution when 

approximations ([5]) and m are reached, we have Q ~ 2M'^ec or 2M'^ee and ~ 1. However, for the region near the 
bounce, evolution are very complicated, and numerical calculations are required. Here we numerically plot Q/c^ and 
Cg in Fig. m 
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FIG. 4: Plots of functions Q/(?s and during the whole evolution. The parameters are chosen as before. We can see that 
during the whole evolution Q/c^ > 0, indicating that we don’t have any ghost instability in our model. In regions away from 
the bounce, Q/c?s reduces to constant proportional to the slow-roll parameter, is trivially equal to unity in regions away 
from the bounce, however around bouncing point, it goes down to below zero, causing short period of gradient instability. 


From the Fig. |4]we can see that, in the whole evolution (including the bounce) we have positive Q/c^, indicating 
that our model is completely ghost-free. However, during the bounce region, the numerics show that there is a instant 
period of negative c^, which may induce fast increase of the perturbations. We will later show what sequent it will 
cause to have negative and how it can be ameliorated. In fact, the same issue also appears in the model given by 
[30l|. Let’s first write down the equation of motion for C, according to lagrangian (HU): 


where 


u" + c^V^u - u = 0 , 


_ _ a^/Q 

u = zQ ^ z = - 


(19) 


( 20 ) 


and “/” and “V” denotes derivatives with respect to conformal time r] = J a~^{t)dt and space x respectively. When 
transformed into momentum space, one has u{x,t) —>■ Uk{t) for each wavenumber k, and can be substituted into 
— k^. According to the numerical result in the last section, we can plot the evolutions of fluctuation modes for different 
wavelengths Xphys = a/k as well as the Hubble horizon of our model as in Fig. [5l From the sketch plot one can 
see that, these fluctuation modes can been divided into two parts, one is of small fc’s (green, cyan and blue), while 
the other is of large fc’s (navy and purple). The former modes are generated in the contracting phase, and will exit 
the horizon before bouncing, while the latter ones are generated in expanding phase, and will exit the horizon after 
bouncing. As will be shown very soon, they will have very different evolution and fc-dependence. We remark the 
dividing mode between the two as fcg (the blue one), the wavelength of which is just equal to horizon at the beginning 
of the bounce region. 
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FIG. 5: The sketch plot of the horizon of our model (red) and the fluctuation modes that cross the horizon before or during 
inflation period (from green to purple). Since before the bounce the universe evolves slowly with a small value of Hubble 
parameter, the Hubble horizon is large at most of the time, with a sharp damp near the bounce. However, when bounce 
happens, it goes to infinity as H = 0. At the inflation era when H becomes nearly a constant, the horizon also behaves as 
a constant. The fluctuations with k < ko exit the horizon at contracting phase, thus will stay outside the horizon all over 
the inflation region. Those with k > ko exits the horizon at inflationary era, and will stay inside the horizon till the universe 
expands, and exit the horizon at somewhere during the inflation era. 


1. initial condition 


Since during the whole evolution, the background is different at each stage, and will lead to different evolution of 
perturbation. In the following subsections, we will analyze evolution in each stage separately, and joint them together 
using matching conditions to get the final power spectrum. 

At far past, every mode is inside the horizon. In this region, we solve the Eq. (IT^ and taking the >> z"/z limit, 
and obtain the solution; 

= . ( 21 ) 


As we know that in contracting phase except near the bouncing phase, Q and Cg can be approximated as 1. Moreover, 
the coefficient A can be normalized by requiring a quantum origin of the perturbations [s^. We assume that the 
perturbations starts from Bunch-Davies vacuum as usual, and from the quantum condition, we have A = (2fc)“^/^. 
So the initial condition for u is: 


Uk 


_}_ ik-q 


( 22 ) 


2. contracting phase 


In the most part of contracting phase, the model will be governed by the potential term while the higher derivative 
terms are severely suppressed. In this case, the approximation ([5]) applies, and from (IT71) and (IT^ one has Q ~ 2Mp€c 
and Cg ~ 1. Moreover, if the universe evolves according to the scaling solution of Eq. (nnD with the EoS being as m, 
the scale factor of the universe in contracting phase can be parameterized as 

- (t? - , (23) 

where fjB- = tjb- — [(cc — 1)'Hb-]~^ , and tjb- and Hb- are the conformal time and conformal Hubble parameter at 
the end point of the contracting phase, respectively. The perturbation equation (fTOl) then becomes 


u'l, -\- k'^u 


2 — Ec u 

(ec- 1)2 [rj-fjB-Y 


= 0 . 


(24) 


The solution of the above equation is Hankal function; 


Uk = Vlv - fiB-\[ci{k)J^_{k\r] - 77b-I) + C 2 (fc) J_^_(fc|7? - iys-D] , (25) 

with the coefficients ci(fc) and C 2 {k) to be determined by the initial conditions ([2^ . Here 

3 - Ec 

" 2(Ee - 1) • 


(26) 
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In the region far away from the bounce where all the fluctuation modes are deep inside the horizon, one have 
\z"lz\^ and the solution (1^ can be approximated as: 


Mfc 


Cfe 


~ Ci(fc) 

Cl (A:) 


cos(fc|?? - rys-l) , C 2 ik)\ — sm{k\ri - fjB-\) , 
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V - VB- 

as- V ttCcA; V^b- - j?B- 


cos(fc|77 - 77b-I) , 


C 2 (fc) 

a-B- 



y-VB- \ 

Vb- - f]B-) 


sin(fc|77 - 77b-I) , (27) 


In the region near the bounce, the horizon shrinks rapidly to a minimum value. In this case, some fluctuation modes 
with small fc, namely <C \z'' / z\ will exit the horizon and become superhorizon modes in the contracting phase. The 
approximate behavior of these modes according to the solution (l25ll will be: 


Uk 


Ck ^ 


ci{k) 


2'^-r(z/_ + l)y/k 


(fc|77 - i^B-l)'--! , 


1 C 2 (fc) 


2-'^-r(-B_ + l)Vk 


(fc|77 - i7B-|)'--i , 


2'^-r(z/_ + l)^/2ecClB— 




2 ^ r(— jy— + l)y/2ec(iB— 


while the large k mode still remains the same as (EH). 

The subhorizon solution in the contracting phase, which describes the modes deep inside the horizon, can be 
connected to the initial condition (|2^ to find out the coefficients ci(fc) and C 2 (fc). Using (1221) and (l27l) we will get: 


ci(*) = ^ , C2(fc) = . (29) 

One can see that while subhorizon modes shows oscillating behavior, the superhorizon modes splits into two 
branches, one keeps constant and the other is decaying. Assuming the two modes are generated equally, one could 
easily understood that for modes that exits horizon at contracting phase, the constant mode will dominate over 
the decaying one, and will be inherited by the perturbations through the bounce and in expanding phase. We will 
investigate the following evolution of perturbations in the coming sections. 


3. through the bounce 


After the Hubble parameter of the universe reaches the minimum value, it will no longer evolve as an monotonic 
function which keeps negative, but rather turns around towards positive values. This will cause the universe contract 
more and more slowly, while finally it stops contracting and becomes expanding, triggering the bounce. The Hubble 
parameter will go from its minimum value to its maximum value, the process during which we call as a bouncing 
phase. During the bouncing phase, the Hubble parameter is increasing rather than decreasing, and shape functions 
such as k{(j)), t{(j)) and G{X, </>), as well as Q and c^, will have nontrivial features, which will make Eq. (IT^ hard to 
solve. For sake of analyticity, we assume that the Hubble parameter increases as: 

H ~ a{t — Ib) , (30) 


when ^B is the time point when the bounce happens. Then we have 


a ~ aBe“^‘ ~ - tBf‘\ 


(31) 


for the scale factor. Here ob refers to the scale factor at Ib- From the numerical calculation, it can be identified that 
around the bouncing point the first two terms of the Lagrangian 0 will be dominant, so one roughly have: 




2M^ki(j))X 
{M^H - GxX^y 




k{(l>) 

6t((f>)X 


(32) 


respectively. Using Friedmann equations, one can furtherly simplify to be 


c ~ - 


11 + 2a{t — tsY 


21+ ?,a{t-tBY ’ 


(33) 



















9 


and near the bouncing point where the t-dependent terms can be viewed as higher order. Moreover, near the bouncing 
point, the bubble parameter will be small compared to the term GxXcj) in the dominator of Q (at bouncing point 
H = 0), where Gx is obtained from G{X,(j)) in dH). From ([6|), one can notice that in our model Gx = -^p(7/fo)i(<(>), 
so Q can be expressed as 

2M4fc(0)X ulX 

^ ~ G^X^r ~ ~ 1^0, ■ 

It is convenient to approximately parameterize X as function of t : 

X ~ , (35) 


where x is the parameter and is the value of X at the bouncing point. From the numerical calculation we find 
that the parametrization can efficiently mimic the evolution of X near the bouncing point. 

It is convenient to transfer everything into conformal coordinate where the time is denoted by the conformal time 
r]. From the definition of rj one can see that for leading order rj — tje = cLg^it — ts), where rys is the conformal time 
at the bouncing point. Then from Eq. (1321) one can express Q and w.r.t. rj approximately as: 


Q 


^'^qXb ^-^al(n-nB? 
■j^a 


and the parameter z turns out to be: 


11 + 2aa%{r] - tibY 
2 1 + 3aa%{r] - tibY 


(36) 


z = a 


^/Q 


— aB\ 


I mix 


"f^a 




(37) 


Though Q/Y > 0 insures there is no ghost instability, < 0 leads to a gradient instability, which mainly affects the 
small-scale perturbation modes. These small-scale modes are within the horizon, i.e. the wavelength A < 'Hb-,'Hb+ 
(noting that the superhorizon modes are hardly affected, see [s^.), in which case the perturbation modes present 
quantum-behavior. At classical level, it still remains an open issue of how to estimate the effect of this gradient 
instability. Here, we will not involve this issue in this paper, which is beyond our scope. 

However, the nasty instability might be avoided if we expect some additional corrections in the high energy region 
such as around the bounce [d^. For example, we phenomenologically assume that there is an additional term in the 
action of gravity such as: 

= i y dri(Yxa^^{t)R^Y^ , (38) 


(3) 

the physics of which is yet to be known, where the 3-metric Ri ^ can be viewed as the scalar component of potential 
term of the 4D gravity. In order not to affect gravity at other regions of evolution, ^(f) is required to go close to zero 
at past or future limit but reach its extreme value at the bounce. Since as a spatial curvature, r!'P ^ (^C)^ only, so 
this term will do nothing to the background evolution as well as the time derivative of Cj but will add another spatial 
derivative term of C such as 


1 

2 


dijcfixa^ ( Q^{t) 


2Mim 


M'^H-GxXcj)^ 


(dcr 


so adding this term, the effective sound speed squared becomes: 


(39) 


c. = c. 


I + ^{t) + 




Q{M^H - GxXY) 


In order to make > 0 around the bounce, we only need 


(40) 


1 + ^(t) 


2M^at) 


Q{M^H - GxXct>) 


< 0 


(41) 


which can be easily satisfied by choosing proper form of ^{t). As an example, we parameterize ^{t) as ^{t) = , 

and plot the evolution of cl(t) in Fig. [51 From the plot we can see that now the sound speed squared will be positive 
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FIG. 6: Plot of the new sound speed squared under correction (I38II during bouncing phase. With the correction, the sound 
speed squared can be kept above zero, avoiding the gradient instability. Various examples are shown with set to be — 1 
(blue), —1.5 (pink) and —2 (gold). At the edge of bounce region where the correction term starts to decay, it will be matched 
to values of in contracting/expanding region, which is unity. 


all the time, avoiding the gradient instability around the bounce, while reduce to the one in contracting phase due to 
the quick decay of ^ (t) far from the bounce. 

The EoM of perturbation now becomes 


u'k + -{a- x)a%) Uk = 0 , (42) 

with a positive c^. As explicit examples, we set = ~1) —1.5 and —2 while Aj = 1, and find that can be easily 

kept above zero (the minimum required value of is —1, which gives Cg ~ 0). The solution is: 

Uk = diCOsll(r/- t/b)] + d2sin[l(i] - t/b)] (43) 

for any k, where P = c^k'^ + {a — x)ag. In general, we have x > a during bouncing phase, and this is the requirement 
of kinetic dominance in order to get a bounce: (j) > H(j). This relation can also be numerically checked in our model. 

The solution during bouncing phase should be matched onto the solution in contrac ting phase via the matching 
conditions. We apply the Hwang-Vishniac (Deruelle-Mukhanov) matching conditions [3^1^ that Cfe and its derivative 
should be matched at the pivot point rjB- (see also [4l[ for matching conditions in Galileon theories.). Using the 
solutions (EZl), dMl) and (031) we have: 


di = I sm[csk{riB - Vb-)] (c 2 cos[fc(?7B- - Vb-)] - ci sm[k{T]B- - j?B-)]) 

+Cs cos[csk{'qB - VB-)] (ci cos[A:(r7B_ - rys-)] -b C 2 sin[A;(77B_ - iys-)]) | , 
d2 = y^ ^g^^2fc | cos[csA:(r7B - Vb-)] (c2 cos[fc(ryB_ - rys-)] - Ci sin[k{r]B- - Vb-)]^ 
-Cs sm[csk{r]B - Vb-)] (ci cos[fc(?7B- - Vb-)] + C 2 siii[k{r]B- - ??b-)]) > 


for subhorizon modes while with (|28)) . (l29l) and (|43l) we have: 

[(ec- 1 )'Hb-]^ 


di = 


2 '^ r(i + B_)y^ 


k 2 + cos[l{-r]B + r]B-)]ci , 


r(i^+1")^ 


(44) 


(45) 


for superhorizon modes. 
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4- expanding phase 


After the bounce, the universe will enter into an inflationary expanding phase, with the Hubble parameter being 
a little bit decreasing but nearly steady. In this phase, all the features of the shape functions will decay again and 
make them the same constant value as before the bounce. Moreover, the parameters Q and will become constant 
again, with Q ~ 2Mpee and ~ 1. The scale factor during this time can be parameterized as: 

aiv) - iv - VB+)^ , (46) 

while fjB+ = tib+ — [(ce — to protect the continuity of "H, and the slow-roll parameter Ce is defined in Eq. (1141) . 

However, since during the most of the time (j) is much smaller than the expectation value v, one have Ce <C 1 and 
ie Hee, thus it is reasonable to take Ce as nearly a constant. Then the equation of motion (11911 becomes: 

u'l + k'^u-j ^—^—12 = > ( 4 '^) 

{I - ie? [r] - Vb+Y 

and the solution is the Hankal function: 

Uk = - - Vb+\) + 92{k)J-u+{k\ri - Vb+\)] . (48) 

where 

3 - Er 


= 


2 (ee - 1 ) ■ 


(49) 


As has been discussed in the subsection 2, the solution of Eq. (I48p also has two branches of solution. One is the 
subhorizon solution: 


Uk ^ 5i(fc) 
9i{k) 


cos{k\r] - fjB+\) : 92ik)\ —sm{k\ri-fiB+\) 


Cfe - 


a_B+ V TTC. 


V - VB+ 
VB+ - ?7B+ 


nk 

fu ~ l^ 92ik) 
cos{k\T] - riB+\) , 


UB+ V Tre, 


'n-VB+ 

r]B+ - fjB+ 


sin(fc|r 7 - 773 +1) ■ (50) 


This solution is for the large k modes, which still remains inside the horizon at least at the time when the universe 
enters the expanding phase. This branch should be matched with the large k modes in bouncing phase (j43ll . The 
other branch of the solution of (H51) is the superhorizon solution: 


Uk 


Cfc 


9iik) 


2^+r(i/+ + l)Vk 

gi(fc)[(ee - 1 )'Hb+]^^ 
2‘"+r(jz+ + l)^aB+ 


(fc|?7 - rys+D'e-i , 


92{k) 


_ €e —3 

-k 2(ee-l) 


2-‘^+r(-i/+ + i)Vk 

92{k)[{ee - 1)Hb+]^ 


2 + l)y/^aB+ 


ik\v-flB+\)‘‘-^ , 

k~ (^k\ri - fiB+\)"'‘^~^ , (51) 


This solution has two origins: one of the origin is from the large k modes, which exit the horizon during the expanding 
phase. Since it shares the same equation with the subhorizon solution (1501) . the coefficients of the two parts will be 
the same. The other origin is from the small k modes, which already exit the horizon when the universe enters the 
expanding phase. This solution should be matched with the small k modes in bouncing phase (1441) . So we have: 


gi = 


1 


92 = 


|2cs cos{cskAr]B)(ci cos[k{r]B- - fjB- - 77 B+ + j?B+)] + C2 sin[fc(?7B_ - fjB- - 7?b+ + ?7B+)]) 
-(Cg + 1) sin(csfcA77B) {c 2 cos[k{ir]B- - Vb- - riB+ + 77 B+)] - ci sin[A:(77B_ - 77 B- - ?7b+ + j7B+)]) 

+(c^ - 1) sin(csfcA77B) (c 2 cos[fc(?7B- - Vb- + riB+ - f?B+)] - ci sin[A:(77B- - Vb- + r]B+ - riB+)]j | , 

2csCOs{cskAr]B)(c2 cos[k{r]B- - VB+ - Vb- + j?B+)] - ci sin[fc(?7B- - ??b+ - Vb- + Vb+ 
2Cs V Cc I ^ 

+(Cg + 1) sin(csfcAr7B) (ci cos[k{riB- - 'nB+ - Vb- + j?B+)] + C 2 sin[A:(77B- - 7 ?b+ - Vb- + 7?b+)]) 

+(Cg - 1) sin(csfcAr7B) (ci cos[k{riB- + r]B+ - f)B- - j?B+)] + C 2 sin[A:(77B- + 7 ?b+ - j?B- - f?B+)]) | 


(52) 
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for those which exit the horizon in expanding phase, where we define Atjb = 'r]B+ — Vb-, and 


9i = 


92 = 


(ee - 3)r(l + u_) V [(e^ _ l)nB+]^ 

- iy+) 

(ee - 3)r(l + 


^[(Ce - 3 )cos(;A?7b) +;'HB+sin(;A?7B)]ciA:'c-i + i-=e , 


— [(Cc - [(Ce - l)'HB+]i-'=Z(ee - 1) sin(;A77B)cifc 


-i+, 


(53) 


for those which exit the horizon before the bounce. 

We are observing superhorizon modes, namely. For the solution of the superhorizon mode (1511) which we’re observ¬ 
ing, the constant branch gi will be dominant over the decaying branch 52 - The power spectrum is then 


^ 27r2'^' 27r2' 2 '^+r(j^+ -b l)V2iIaB+ ' 


(54) 


Considering that ci = ci = iy/TTj2, the modulus squared of gi{k) is calculated as: 


l5i(fc)P = 


TT £e 
16c2 Cc 


-( 1 - 


-2 


(1 -b Cg) -b(l —c^)cos 
2k 


2k 


- (1 - ci) cos{2cskAr]B) 


(1 -b ci) cos 


(Ce — 1)'Hb+ 


(€e - mB+ 

cos(2cskAgB) + 2cs sin 


2k 


(€e - mB+ 


sm(2cskAgB) 


>(55) 


For large k modes, the trigonometric functions in the above expression will oscillate rapidly, showing only their average 
values, which is zero. Then the mean value of \gi{k)\^ ~ (7r/16Cg)(ee/ec)(l -b c^)^. 

One can see that, for modes that exit horizon before bounce. 


(cc-l)i-'-H|_[(ee-3)cos(;A77B)-b;H5+sin(lA7?B)]2 / k \^cA 
22--+47r(ee-3)2r2(l + j2_)e, \n^) 


with the index 


SCf- 


- 1 


Since Cc > 3, it will have a blue tilt. 

While for modes that exit horizon after bounce, the power spectrum F’^(fc) is in the form of: 

^ (l+c1)^(ee-l)T^g|^ / k \ ^ 

^ 647recc222-+r2(j2+ + l) \nB+ J 


and the spectral index 


n 


inf 

C 



(56) 


(57) 


(58) 


(59) 


We can see that, taking the correction (l38l) . it shows scale-invariance even at the fc —>■ 00 limit, which is consistent 
with the data. When we take the slow-roll limit Ce — 0, we have ~ (1 -b Cg)^i7^_|_/(327r^ecCg). Moreover, if one also 
take into account the corrections from the trigonometric functions in \gi{k)\'^, both P^ and will have some wiggling 
behavior due to such a correction, which will be observed if we have enough accurate data. The same property has 
been discussed in the second paper of Ref. m- This property can be used to distinguish from the standard inflation 
models without preceding evolution. 


B. tensor perturbations 

Besides the scalar part, there are also tensor part perturbations. Since the field is pure scalar field and there are 
no tensor degrees of freedom in the field part, the tensor part totally come from the perturbations of the gravity. 
It is believed that the tensor perturbation is originated from the primordial gravitational waves, and although very 
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weak, the recent gravitational wave explorers such as BICEP has shed some light on our future ability of detecting 
them more clearly . As a theoretical study, here we also discuss the tensor perturbations generated in our model. 
According to [s^, the equation of motion for tensor perturbations can be written as: 


c( 2 ) _ 1 

- 2 


drjd^xa^ 


h'^ - {dhf 


(60) 


where h is the amplitude of the spatial-spatial part perturbations of the metric, namely hij. There are two degrees of 
freedom in tensor parts, namely and obeying the same action and equation of motion. From (1601) we get the 
equation of motion for h as: 


< + (fc" --W = 0, 


(61) 


where we define v = ahj2. From both the action and the equation of motion we can see that, the evolution of the 
tensor part will be affected only by the evolution of the universe itself, namely the scale factor a, rather than the 
functions containing (j). This is due to that tensor part is decoupled from scalar part, and only comes from the gravity. 
Thus solving the equation of motion for tensor perturbations will be much easier than that for scalar ones. 

We follow the same procedure as what we did in subsec. A. At the very beginning when all the modes are deep 
inside the horizon, the last term in the left hand side of Eq. dSB) is not important, thus similar to scalar part, v 
has the oscillating solution. For the initial condition of h, we also choose the Bunch-Davies vacuum, so the tensor 
perturbation at the initial time have the solution as: 


1 


Vk = 




Jkrj 


(62) 


For contracting phase where we have 0 ( 77 ) ^ {rj — fjB-) , we have the solution for tensor perturbations as: 

Vk = \/\'n- fiB-\[cl{k)j^_ ik\r] - fjB-\) + cl{k)J_^_ {k\r^ - 7 ) 5 -1)] , (63) 


where v- is already given in Eq. (1261) . Since in this period the functions such as Q and for scalar perturbations also 
becomes trivial, the solution of tensor and scalar perturbations are actually very similar. For large k modes which 
remains inside horizon during contraction phase, the solution is: 


Vk c((A:)W —cos(fc| 77 - 7 )B_|) , C 2 (fc)W —sin(fc |77 - 77 B_|) , 


hk — 


cT(k) 

Qb- 


1 - VB- 
r]B- - Vb- 


cos(fc|77 - fjB-\) , 


Kk) 


a-B- 


V-VB- 


riB- 


riB- 


sm{k\r] - fjB-l) , (64) 


while for small k modes which has already exit the horizon during contracting phase, the solution will be 


Vk ^ 


hk ^ 


zj{k) 


2^-r{v_ + l)Vk 

cj{k)[{ec - 1 )'Hb 


{k\r] - i 7 B_|)'c-i 


\k) 


-k 2 ('c-l) 


2-''-r{-v_ -f l)^/k 
C2’(fc)[(ec - 1)'Hb-]~ 


2''--^r{i^_ + l)aB 
By connecting the subhorizon solution with initial condition 


(fc|77 - i7B_|)-c-i , 

_ €0 —3 ec —3 

k 2(ee-l) (^k\r] — 77 B- I) • 


- 
Cl — 


ir(-i/_ -I- l)aB 

, one can get the coefficients: 


T 

C 2 = I 


(65) 


( 66 ) 


For bouncing phase however, the tensor perturbation shows greatest difference with the scalar perturbations, for 
it is independent of functions such as Q and c^, which becomes nontrivial any more in bouncing phase. Considering 
Eq. (EH and the expression of scale factor in bouncing phase (EH, the equation of motion for Vk now becomes (in its 
leading order): 


and the solution is: 


v'k + (fc^ - aa%)vk = 0 , 

Vk = cos[;(77 - t^b)] + sin[Z(?7 - 77 b)] , k > ^JaaB , 

Vk = , k < y/^aB , 


(67) 

( 68 ) 
(69) 
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where P = \k‘^ — aa^|. Note that rj will run from rjB- to rjB-, and the value of Ufc at 77 = rjB- is the amplitude growth 
of the tensor perturbations during bouncing period. Via matching conditions which requires continuity of hk and 
l4nl |. one can determine the coefficient as: 


dn 




[c\ cos{k{r]B - VB-)) + C 2 sm{k{r]B - Vb-))] 


[c 2 cos{k{r]B - Vb-)) - cf sin(A:(r7_B - Vb-))] 


(70) 


for subhorizon modes, while 


^ cn(ec - mB-] 

i Ol-l-7> 7T'/1 I .. \ ' 




2^+''-lT{l + v_) 


(71) 


for superhorizon modes. 

In expanding phase, the tensor perturbations becomes concordance with the scalar ones again. Since the scale 
factor now becomes 0 ( 77 ) ^ (77 — fiB+Y^^^~'^‘\ the solution for Vk becomes: 


'^k = \/\{v - TlB+)\[ 9 l[k)J^^{k\ri - 77B+I) +5^ {k)J_^^{k\ri - 773 +!)] , 

Also the solutions are divided by subhorizon and superhorizon parts. The subhorizon part is: 


(72) 


Vk ^ (A:)W — cos(fc|77-?7 b+|) , 52 (fc)v^sin(A:|77 - 775+1) 


hk — 


gfik) f 77 - 7 ) 5 + V"'" 


05 + V nk \r]B+ — fiB+ 


Tvk 

fU ~ |N 92 ik) [Y ( T]-fjB+ 

COs(fc|77 - 775+|) , - \ —r - 1 - 

05 + V nk \r]B+ - r]B+ J 


sm{k\r] - fjB+\) . (73) 


which denotes the modes that remains inside the horizon at the beginning of the expanding phase. While the 
superhorizon parts is: 


Vk ^ 


aTik) ,,, . |,_ L _^ 5j(*) 

——- 77 - 775 +)'<=-! , ;7 -—^(A: 77 - 775+ , 

2''+T{v++ l)Vk 2-''+T{-v+ + l)Vk 


+ ^ , gSy , (T4) 


2’'+ W( 7 z+ + 1 ) 05 + 


2 ''+ ir(— 17 + + 1 ) 05 + 


which denotes the modes that are outside the horizon. These modes also have two origins. For those which exit the 
horizon in expanding phase, using matching conditions with the solution in (16811 and (I73|) . we have: 


(75) 


gj ~ Cl cos(A:(775+ - fjB-)) + C2 sin(fc(775+ - fjB-)) , 
g^ ~ C2 cos(fc(775+ - fjB-)) - Cl sin(fc(775+ - fjB-)) , 


and for those which exit the horizon before the bounce, one obtain 


2"+-‘^-r(l + 77 +) [(e, - 1)775-] 

9i - 


T 

92 - 


{ce - 3)r(i + 77_) _ i)nB+]^‘ 

2-''+-''-r{i - 77+) 


(Ce - 3) cosh(ZA775) + —-sinh(ZA775) ) Cifc'^-i + i-'-^ , 

rLB+ 


(ee - 3)r(l + 77_) 


[(Cc - 1 ) 775 -] [(Ce - 1 ) 775 +] l{^e - 1) sinh(;A775)cifc 


- 1 +- 


(76) 


Similar as the scalar part, we are observing superhorizon modes, namely dZH). For this solution, it has two branches, 
one is constant ( 51 ) and the other is decaying 52 , therefore the first one is dominant. The power spectrum is then 




= k\ 9l{k)[{e,-l)nB+V-‘ 

2'"+r(77+ + 1 ) 05 + 


( 77 ) 
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where the factor 2 appears because there are 2 independent polarizations of the graviton [i^. One can see that, for 
modes that exit horizon in the inflation phase, 


(ce - 1) / k \ 

22-+^r2(z.+ + l) \Hb+J 


(78) 


with the index tit — 2ee/(ee ~ !)• For Cg — 0, we have Pt = 2{Hb+/tt)‘^ and nr — 0. While for modes that exit 
horizon before bounce, 


Pt = [(ec-1)]^ 


2 2 [(ce - 3) cosh(lAryB) + sinh(/A?7B)]2 / k 




22^-7r(eg-3)2r2(l + i/_) 


\'Hb- 


(79) 


with the index nx — 2ecl{ec — !)• Since Cg > 3, it will have a blue tilt. 

From the above result for scalar and tensor perturbations, we can also calculate the tensor-to-scalar ratio r, defined 
as r = Pt/P(;. The 1st year Planck paper suggest to have r < 0.1 at 2cr level, while BICEP2 data shed some light on 
the possibility to measure it more accurately. In our model, we have 


64ecC^ 

(1+C2)2 


for large k , 


r = IGcg for small k . 


(80) 


One can see that for large k modes which is observable, r carries information from bouncing or contracting phase, 
which makes bounce detectable. For example, the Planck constraint r < 0.1 requires that < 5.2 x 10“"^ for the 
minimum value of eg = 3, which furtherly lead to 0 > (1 + ^o) > ~3 x 10“^ for in the parametrization of ^{t). 
However, for small k modes which has not been detectable yet, r could be very large (> 50). This could be a smoking 
gun to the bounce inflation scenario. 

As a side remark, we would like to comment that actually all the above analysis is based on the consideration that 
the bounce is not far from observable inflation, i.e. the e-fold number is not too much, so the information from the 
bounce can be detectable. However, there is also possibility that if the inflation was too long, therefore the observed 
modes comes from Bunch-Davies vacuum even in expanding time (like the purple line in the Fig. [5]). In this case, 
we instead have gi = \pri 12, g 2 = 1-07/2, so everything could be completely the same as pure inflation and bounce 
becomes non-detectable. As for in the region before the bounce, this mode should be inside quantum phase all the 
time, and solutions in different stages may be matched quantum mechanically. Of course such an issue goes far beyond 
the scope of this paper, and actually we don’t know how to do it yet. There may be chances that we pick up this 
topic in future investigations. 


IV. FITTING THE DATA 


In this section, we try to fit our model with the newest released PLANCK data. Following (1551) and 
parametrize the whole power spectrum as: 


we 


ppara _ 


1 — tanh 




XTLb- 


- 1 


^cou 


T-Lb- 


+ 


1 -I- tanh 


T-Lb- 


- 1 




'Hb+ 


(81) 


where A^°'^ and are the fc-independent prefactors of expressions (1561) and ([58|) . and n/°” and refers to 
Eqs. (l57)) and ((59l) respectively. According to this parametrization, we plot the CMB TT power spectrum of pure 
power-law inflation model and our model in Fig. [71 From the figure we can see that, in most of the plotted region 
{I > 10) our model is identified with the inflation model, and fits the data very well. This is actually not strange, 
because the perturbations of observable scales exits the horizon during inflation era, and should carry information 
of inflation part. The tiny errorbars of the data in this region also indicate that the measurement has been done 
accurate enough, and the shape is almost confirmed and hardly to have other alternatives. However, the data at very 
large scales has large error bars, and the data is not so accurate. Moreover, the data points mildly favor that there is 
suppression in this region. In pure power-law inflation models where all the fluctuation modes are scale-invariant, it 
is difficult to obtain such an suppression, however, in our models where the fluctuation modes exit the horizon before 
bounce, suppression takes place due to the blue-tilt of the spectrum. This result is consistent with [H, [2^ , which 
could be seen as an important distinct from our model and inflation. 
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FIG. 7: Best-fit TT power spectra for pure power law (green solid) and our model (dashed) using Planck-|-WP data. The 
parameter choices are: = 0.02203, = 0.1204, r = 0.09, 1000s = 1.04, As = 2.1 x 10“®, ris = 0.961, = 0.02. The 

red points show the Planck data with la errors. 


V. CONCLUSION 

Inflation theories, although has achieved great accomplishments, has actually been suffering from the singularity 
problem. Phenomenologically, considering a non-singular bounce happens before inflation can help avoid the singu¬ 
larity. The bounce indicates that prior to our expanding universe, the universe might have been contracting from a 
large volume. The bounce scenario not only provide an answer to the mysterious question of “what is our Universe 
like even before inflation”, but more scientifically, may bring some new features to the early stage of our Universe, 
and provide predictions that can be discovered in the future observations. 

However, to add a bounce before inflation, one should care about whether it could bring new problems. To our 
knowledge, the ghost instability and the cosmic anisotropy problems are among the most serious problems of the 
bounce scenario. One is due to the NEC violation which is necessary for the bounce, and the other is due to the 
growth of anisotropy during contraction. In this paper, we present a model making use of the Galileon theory, which 
can violate the NEC and realize bounce behavior without introducing ghosts, because of the delicate design of the 
Lagrangian. Moreover, a negative potential in the contracting phase, which causes the EoS of this model large than 
unity, is useful to avoid the dominance of the cosmic anisotropy. In order to enter into the usual inflationary era 
after the bounce, the potential of the model will have to be connected to some positive value functions, which can be 
realized phenomenologically by shape functions. 

Since now there are multi-stages in the evolution of the early universe, primordial perturbations can be generated in 
either in contracting phase or expanding phase. Those generated in expanding phase corresponds to large wavenumber 
k (small scale), which has entered into horizon and can be observed by nowadays observations. Taking proper inflation 
potential, it is not difficult to get nearly scale-invariant Us with proper tilt, which is consistent with the PLANCK 
data. Meanwhile, those generated in contracting phase corresponds to small k (large scale), which reenter into horizon 
much later, and can only be observed by future experiments. Due to the large background EoS in the contracting 
phase, the spectrum of these modes are blue-tilted, which will be more and more suppressed when the scales get larger 
and larger. However, there are already some hints in nowadays observation, i.e., the small I suppression observed at 
TT spectra of CMB map. Moreover, the spectrum might present some oscillating features due to the bounce effect, 
and this can be viewed as a distinction to pure inflation models. Tensor perturbations and tensor spectrum are also 
discussed in this paper, which has similar behaviors as the scalar ones. 

Nowadays along with the rapid development of the astrophysical observations, the observable range are getting 
more and more enlarged, and the observational data are also more and more accurate. More and more interesting 
new phenomena are discovered, which is not only chance but also challenge to the early universe model buildings. For 
instance, the newly released PLANCK data has pointed out that there have been several anomalies in the CMB sky, 
which has been stimulating great interests among astrophysicists and cosmologists. Can the bounce inflation scenario 
be good enough to explain those phenomena? Or what more phenomena can it predict that could be discovered by 
the future observations? We are trying to answer these questions in our upcoming works. 
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